Astron. Nachr. / AN 999, No. 88, 789 -17951(201 1) / DOI please set DOI! 



Transport of angular momentum and chemical species by anisotropic 
mixing in stellar radiative interiors 

L.L. Kitchatinov 1 2 3 * and A. Brandenburg 3 4 

1 Institute for Solar-Terrestrial Physics, P.O. Box 291, Irkutsk 664033, Russia 

2 Pulkovo Astronomical Observatory, St. Petersburg 196140, Russia 

3 Nordita, AlbaNova University Center, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden 

4 Department of Astronomy, Stockholm University, SE-10691 Stockholm, Sweden 

The dates of receipt and acceptance should be inserted later 

Key words stars: interiors - stars: rotation - stars: abundances - hydrodynamics - turbulence 

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., instability of rotational shear. In this 
paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification 
induces strong anisotropy with a very small ratio of radial-to-horizontal turbulence intensities. Angular momentum is 
transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. 
This non-diffusive transport known as the A-effect has outward direction in radius and is much more efficient compared to 
the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed 
using direct numerical simulations combined with the test-scalar method. As a consequence of the non-diffusive transport 
of angular momentum, the estimated characteristic time of rotational coupling ( < 100 Myr) between radiative core and 
convective envelope in young solar-type stars is much shorter compared to the time-scale of Lithium depletion (~ 1 Gyr). 
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1 Introduction 

Mixing of material in stably stratified stellar radiative inte- 
riors is opposed by buoyancy forces. Nevertheless, a certain 
level of turbulence resulting, e.g., from instabilities of rota- 
tional shear (Goldreich & Schubert 1 19671 1 or forcing from 
the upper convection zone (Blocker et al. [1998 ) is possible 
and even necessary to explain rotational coupling between 
the radiative core and the convective envelope, and thus the 
depletion of light elements in solar- type stars. 

Turbulent transport in stellar radiation zones is supposed 
to be strongly anisotropic (Spiegel & Zahn 119921 Denis- 
senkov 1201 Oi l. Buoyancy forces suppress radial mixing so 
that much higher horizontal compared with radial velocities 
might be expected. The degree of the anisotropy is not cer- 
tain, however. It remains a free parameter of transport mod- 
els. The anisotropy in rotating fluids is, however, not free 
(Rudiger & Pipin BOOTT l. The influence of the Coriolis force 
on horizontal motions produces radial displacements. The 
ratio of radial to azimuthal mixing intensities is controlled 
by the balance between Coriolis and buoyancy forces. 

The aim of this paper is to estimate the effects of rota- 
tion and stable stratification on turbulence. We do not spec- 
ify here the origin of the turbulence but just prescribe the 
so-called 'original turbulence' that would take place in non- 
rotating and neutrally stratified fluid. The combined effect 
of rotation and stable stratification on this original turbu- 
lence is then estimated. 
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We shall see that the resulting parameters of turbulent 
mixing depend only slightly on the properties of the pre- 
scribed original turbulence. In particular, the anisotropy, 

(u 2 r ) ~ r> 2 ' K ' 

varies only slightly between very different prescriptions (in 
this equation, u is the fluctuating velocity, fl is the angu- 
lar velocity, N is the buoyancy frequency, and r is the eddy 
turnover time). Typically, £1 <C N in radiation zones and 
probably we also have tN 3> 1. Therefore, the anisotropy 
of Eq. (03 is high. Our main finding is that the strongly 
anisotropic turbulence is much more efficient in transport- 
ing angular momentum than chemical species. The reason 
is that the Coriolis force makes the azimuthal and radial 
motions correlated. Fluid particles with positive (negative) 
azimuthal velocity are pushed by the Coriolis force in the 
positive (negative) radial direction. As a result, (urU^) > 0, 
and the turbulence transports angular momentum in the pos- 
itive radial direction. This angular momentum flux cannot 
be interpreted as an effect of eddy viscosity, i.e., it is not 
diffusive by nature. Transport of chemical species, on the 
contrary, is only due to radial turbulent mixing. 

The non-diffusive transport of angular momentum that 
we find is essentially the well-known A-effect of differen- 
tial rotation theory (Lebedinskii |1941j Rudiger [l989l >. The 
effect has been well studied for nearly adiabatically strati- 
fied stellar convection zones. Our paper suggests that it can 
be important for the rotational coupling between radiative 
interiors and convective envelopes as well. 

© 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



790 



L.L. Kitchatinov & A. Brandenburg: Turbulent transport in stellar radiative cores 



2 Mathematical formulation 

2.1 Relation between the original and background 
turbulence 

To specify the effects of rotation and stable stratification on 
turbulence, we follow the standard approach of quasilinear 
theory by formulating the linear relation, itj = D^v/p be- 
tween the velocity fields of original turbulence u^°\ which 
would take place in nonrotating neutrally stratified fluid, and 
actual or background turbulence u. The relation tensor Dy 
includes the effects of rotation and stratification. 

We apply a simple version of the r-approximation to 
write the equations for fluctuating velocity and entropy as 
follows, 



2fi x u - 



1 



r 



-VP-g = f. 



- = -(u-V)<*>, 



(2) 



where / is a random force driving the turbulence, r is the 
eddy turnover time, g is gravity and the r-relaxation terms 
replace the nonlinear terms together with time-derivatives, 



On u 
— + {u- V)«- ((u- V)u) -»• -, 
at t 

— - + u-Vs - (u ■ Vs) -4 — . 

at t 



(3) 



Microscopic diffusion is neglected. The r-approximation 
([3]) assumes that the nonlinear interaction of turbulent ed- 
dies leads to their effective dissipation (a turbulent fragmen- 
tation of scales) in a characteristic time r s» Iju, where I 
is the characteristic spatial scale of the eddies. We assume 
that the motion is incompressible, divw = 0. This means 
that vertical displacements are small compared to the den- 
sity scale height and velocities are small compared to the 
speed of sound. 

Fourier transformation is applied, 



u(r) 



e ikr u(k) 6k, 



s'(r) = J e ik - r s'(k) 6k, 



(4) 



to convert the partial differential equations (f2]i into algebraic 
equations. We neglect variations of mean fields on the spa- 
tial scale of the turbulence to obtain a closed equation for 
the fluctuating velocity: 



u l + N* 2 (fi - (iki)(f ■ u) + aCl* EipmkpUm. = £t (0) 



(5) 



ipm ^p 

(cf. Kitchatinov & Rudiger (2008 ) on how the pressure term 
is treated). In Eq. (0, f is the radial unit vector, fc = fc/fc 
is the unit vector along the wave vector, /x = (fe • f) is the 
cosine of the angle between wave vector and radius, and 
a = (fe - S~2)/r2 is cosine of the angle between the wave 
vector and rotation axis. The Coriolis number 



2rfl 



(6) 



and the normalized buoyancy frequency 

T 2 9 d{s) 



N* = tN, N z 



c p dr 



(7) 



parameterize the effects of rotation and stratification. 
Solving Eq. © for u gives the relation tensor, 
1 

ij = \ + n* 2 {\- ^) + cf 2 n* 2 x 



1 + N* Z {1 - % + a{l*e t]p kp 



N* 2 (n- /j,ki)rj + 



((fi - fiki)t 



(8) 



This equation describes the joint influence of rotation and 
stable stratification on the turbulence. 

In the limit of neutral stratification, N* — > 0, Eq. ® 
reduces to the familiar expression (Kitchatinov 1986) 



D 



(9) 



describing the influence of rotation in the r-approximation. 
In the other limit of slow rotation, ft* — > 0, Eq. ([S]) reduces 
to 

2 



D ^ — 5ij 



N* 



l + iV* 2 (l-/i 2 ) 



(n - fiki)fj. 



(10) 



This equation accounts for the effect of stable stratification 
alone. 



2.2 Original turbulence models 

The original turbulence properties can be prescribed by 



specifying the spectral tensor Q 



(0). 



<ii| 0) (fe)i2f (*')> = S(k + k') 



(o), 



(11) 



(Rudiger 1989 ). The relation tensor ([8]) can then be used to 
account for the effects of rotation and stratification, 

o — n n o(°) 



Qij = (uiUj) = / Qij(k) 6k 



(12) 



For the simplest case of isotropic nonhelical turbulence 
the spectral tensor reads 

S(<i)_E(Je), A /. / i (13) 



(14) 



In this equation, E is the spectrum function 



2\(0) _ 



<«') 



E(k)6k, 



A more realistic model is anisotropic turbulence with a 
preferred direction being the radial one. Horizontal veloci- 
ties are expected to be much larger than the vertical ones. 
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For the extreme case of strictly horizontal random motions, 
we have 



Q$>=q{k,n)x 
((1 - p 2 ){5 io - hkj) - (n- fj,ki)(fj - fikj)). (15) 

The correlation lengths for vertical and horizontal directions 
may differ. If, however, the lengths are equal, the g-function 
in Eq. ( TTBT ) does not depend on p: 

3E{k) 



q = 



(16) 



8vrfc 2 ■ 

The case where the correlation length in radius is small 
compared to the horizontal correlation length can be mod- 
eled by the equation 

3E(k±) k x 



Q = 



8nk 3 



(17) 



The spectrum functions E of Eqs. (fT~6b and ( fTTI i are related 
to the turbulence intensity by the same equation ( TBI as in 
the case of the isotropic turbulence of Eq. JT3V 

Anisotropic turbulence that has finite radial velocities 
(but different from the horizontal velocities) can be modeled 
by a linear superpositions of the spectral tensors dT3l > and 
(fTBT l. We shall see that, if the buoyancy frequency is large 
compared to the rotation frequency, N ^> O, and the nor- 
malized buoyancy frequency is large, N* ^> 1, which are 
conditions typical of stellar radiation zones, then the turbu- 
lent transport parameters differ little between the represen- 
tations ( fT3] >, ( fT6] l and ( fTTl l for the original turbulence, i.e., 
the transport characteristics are not sensitive to a particular 
choice of the original turbulence model. 

2.3 Direct numerical simulations 

An independent verification of the effects of strong stratifi- 
cation on turbulent diffusion is provided by means of direct 
numerical simulations. In that case, we solve the full set of 
compressible hydrodynamic equations for p, it, and s: 

^ = -pV.«, (18) 

P— = -VP+V ■(2p^S)+ P (f + g) 1 (19) 

pT^ = V • KVT + 2pvS 2 - -p(c 2 - c 2 s0 ), (20) 
where D/D£ = d/dt + u ■ V is the advective derivative, 



■|(5yV ■ u is the traceless rate 
of strain tensor, commas indicate partial differentiation, v 
is the kinematic viscosity, the specific entropy is given by 
s = c y In P — c v In p, where c p and c v are the specific heats 
at constant pressure and constant volume, respectively, the 
temperature is related to P and p via (c p — c v )T = P/p, 
which, in turn, is related to the sound speed c s via c 2 
- I ' i). where 7 = c p /c v is the ratio of specific heats, / is 
the external forcing function, and K is the thermal conduc- 
tivity. The last term in the entropy equation (f20b represents 
a cooling term that keeps the temperature (or sound speed) 
approximately constant with a given relaxation or cooling 



time t s . The flow is driven by a random forcing function 
consisting of nonhelical waves with wavenumbers whose 
modulus lie in a narrow band around an average wavenum- 
ber kf (Haugen et al. 120041 1. We arrange the amplitude of 
the forcing function such that the Mach number based on 
the rms velocity remains below 0.1, so the effects of com- 
pressibility are negligible. 

We consider a cubic domain of size L 3 with periodic 
boundary conditions in the x and y directions and in- 
sulating impenetrable stress-free boundary conditions on 
z = ±L/2. The mass in the volume is therefore con- 
served and given by J pd 3 x = poL 3 . The lowest wavenum- 
ber that fits into the domain is k\ — 2w/L. Our average 
forcing kf wavenumber is chosen such that fcf/fci = 5. 
We adopt isothermal initial conditions with c s = c s o and 
p = Po exp(— z/H), where H — c^/jg is the scale height 
which is chosen such that Hk\ = 1. The turnover time 
based on the wavenumber ^ is t = (u Tms kf) , where 
u rms is the rms velocity based on all three velocity com- 
ponents. We vary the normalized buoyancy frequency by 
varying the forcing amplitude. For each run the viscosity is 
adjusted such that the Reynolds number, Re = u ms /vkf, 
is around 40. The Prandtl number, Pr = vc v p/K, is chosen 
to be equal to unity. The cooling time is chosen such that 
TsUrmsfcf is also around unity. 

To quantify the suppression of vertical mixing in a nu- 
merical simulation of stratified turbulence we use the test- 
scalar method by solving equations for the fluctuating pas- 
sive scalar concentration in the presence of a prescribed 
mean passive scalar concentration, 

8r pq 

^— = -V • (uC pq + uc pq - (uc pq )) + n\7 2 c pq , (21) 
at 

where u is the velocity fluctuation obtained from Eq. ( fT9l , 
and the mean flow is zero. The superscript p (= 1 or 3) stands 
for test scalars varying in the x or z directions while q stands 
for c and s that represent mean fields that are proportional 

to cos kx p or sin kx p , respectively, 

C pc = C cos kx p , C ps = C sin kx p , (22) 

where (x\, X2, #3) = (x,y,z) are Cartesian coordinates 
and Co is a normalization factor. Angle brackets denote pla- 
nar averaging over the yz plane (if p = 1) or the xy plane 
(if p = 3). The passive scalar flux is then related to the gra- 
dient of C pc via (uiC pq ) = —Kij'VjC 1 " 2 , where are the 
components of the turbulent diffusivity tensor, obtained as 

k X x = (sin kx{u x c lc ) — coskx(u x c ls ))/kCo, 



k Z z = (sin kz{u z c lc ) — cos kz(u z c ls ))/kCo, 



(23) 
(24) 



by using the solutions for four different test scalars. These 
coefficients depend on time and one spatial coordinate, but 
because of stationarity and approximate homogeneity of the 
turbulence intensity, we present in the following temporal 
and spatial averages of these coefficients. 

The test scalar method has been used previously to 
quantify mixing in turbulence in the presence of rotation and 
magnetic fields (Brandenburg et al. 2009), shear (Madarassy 
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& Brandenburg 2010), as well as isothermal density strati- 
fication (Brandenburg et al. 2012). However, unlike those 
earlier works, the entropy equation is here included and a 
non-isothermal equation of state for a perfect monatomic 
gas is used with 7 = 5/3. Equations (fl~8T> — (f2TT> are solved 
using the PENCIL Codi{] in three dimensions. The numer- 
ical resolution used for the simulations presented here is 
128 3 meshpoints. Like the Prandtl, the Schmidt number, 
Sc = v/k, is also chosen to be equal to unity. 

3 The effect of stable stratification 

We consider first the effect of stable stratification for a non- 
rotating fluid, VL = 0. The only preferred direction in this 
case is the radial one. For any model of the original turbu- 
lence, we therefore have 



(UiUj) 



«) 



(A 



nr.i) + 



(25) 



where (u^) = (u 2 ^) + (u 2 ) is the horizontal turbulence in- 
tensity. Standard spherical coordinates (r, 6, <fi) are used. 

On using the relation tensor ( [Tol l and Eq. (fT2l for the 
isotropic original turbulence of Eq. (1131 . we find 



2\(0) 



AN* 



8N* 3 VN* 



(u 2 r ) = 



In 



+ 1 + N* 



VN* 2 + 1 - N* 



(26) 




- 1 



(27) 



\2N*y/N* 2 + 1 \VN* 2 

Figure Q] shows the ratio of vertical to horizontal turbulence 
intensities as a function of the normalized buoyancy fre- 
quency N*. Squares on the plot show the results of direct 
numerical simulations described in Sect. 12.31 

Horizontal turbulence intensity is influenced by the sta- 
ble stratification only slightly. In the limit of very large N*, 
it is (tt?) = (it 2 )' ) /2, i.e., (it 2 ) is reduced compared to 
the case of a non-stratified fluid by a factor of 3/4 only. The 
vertical motions, by contrast, are strongly suppressed: 

/, ; 2\(0) . . 

(u 2 r ) ~ (ln(4iV* 2 ) - l) for N* » 1 (28) 

in the strong stratification limit. The originally isotropic tur- 
bulence is changed towards horizontal turbulence as N* in- 
creases. 

As might be expected, the horizontal turbulence of Eq. 
( TT3T > is not influenced by the stratification at all, i.e., Qij = 
Q\? in this case. Therefore, whatever model for the original 
turbulence in Section 2.2 is used, the resulting turbulence in 
the limit of large N* is almost the same. The only difference 
between isotropic and horizontal original turbulence models 
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Normalized buoyancy frequency 

Fig. 1 Ratio of the vertical to horizontal turbulence inten- 
sities in dependence on the normalized buoyancy frequency 
(0 for the case of a non-rotating fluid and isotropic original 
turbulence. Squares show the results of 3D direct numerical 
simulations for isotropically forced turbulence. 
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Fig. 2 Dependence of n xx /Kto (open symbols) and 
Kzz/Kto (filled symbols) on the normalized buoyancy fre- 
quency. The dashed line shows that K xx /n t o ~ 2.5 while 
the solid line gives k zz /k w « 0.09 (rA^) _3/2 . 



is that for the first case the radial velocities of Eq. d28l ) are 
finite, although very small. 

Similar results are obtained numerically for the turbu- 
lent diffusivity using the test-scalar method; see Fig. [2] 
where we plot k xx and k zz , both normalized by the ref- 
erence value K t o = «rms/3£;f. Note that K xx /nto ~ 2.5, 
independent of tN, while K zz /n t o « 0.09 (rNy 3/2 . Vi- 
sualizations of the scalar quantity rV • u on the periphery 
of the computational domain (Fig. [3]) confirm our expecta- 
tion that turbulent structures that become flatter as tN is in- 
creased. These structures become smoother in the xy plane, 
while their length scale in the z direction decreases. 



http://www.pencil-code.googlecode.com 
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Fig. 3 Visualizations of rV • it on the periphery of the computational domain for tN — 0.17, 0.9, and 3.7. 



4 Mixing in rotating and stratified fluids 

Similarity of the results for different original turbulence 
models is even more pronounced when rotation is included. 
Coriolis force produces radial motion even for original hor- 
izontal turbulence. The relative intensity of radial mixing is 
controlled by the balance between buoyancy and Coriolis 
forces. 

Producing analytical results for arbitrary values of TV* 
and CI* is problematic. Derivations were performed for the 
practically interesting case of strong stratification, N* ^> 1, 
and not too fast rotation, Cl* N*. The period of gravity 
waves for radiation zones of solar-type stars is of the order 
of one hour. The characteristic time of turbulence is proba- 
bly much longer (it is about one month for convection and 
probably longer for radiation zone motions), so that the con- 
dition N* > 1 is well satisfied. The condition Cl* < N* 
should be satisfied as well in radiation zones of pressure 
supported stars (otherwise we are dealing with centrifugally 
supported disks). 

We now use Eq. ([8]) for the relation tensor that includes 
the combined effects of rotation and stratification. For the 
case of the isotropic original turbulence of Eq. ( fT3l , the 
terms of lowest order in N* 
locity correlations read 
1 

r 

\(°) r / an* 



in the most significant ve- 



«) 



M) 



2\(0) 



AN 



*4 



hi 



l + n*^cos 2 6» 



1 + CI* 



rr 2 cos 2 ' 



9 In 



4N* 



1 + Cl 



' 2 COS 2 I 



(u r u s ) 



= /„ 2 \(°) 



Cl* 



v " ' AN* 4 
x sin 9 cos 9, 
Cl* 



In 



AN* 



l + Cl* 2 cos 2 9 



-3 



= /„ 2 \(°). 



4N* 

The turbulence intensities 



■sin 0. (29) 

z 2 ,) and (it 2 ) control the eddy 
diffusion in latitude and radius, respectively. The cross- 
correlation (urUff,) is important for transport of angular mo- 
mentum. The correlation (u r ue) may cause the temperature 
variation with latitude; positive (u r ug) cos 9 implies pole- 
ward eddy heat flux. It may be noted that rotation does 



not produce anisotropy in the horizontal plane in the strong 
stratification limit, (it 2 ,) = (u 2 ) = (u 2 )/2 (more precisely, 

the anisotropy is small: {(u 2 e ) - {u^))/(uD ~ Cl* 2 /N* 4 ). 

The results for horizontal turbulence of Eq. < fT~5T > with an 
isotropic correlation length of Eq. ( TT6I ) or a short vertical 
correlation length of Eq. ( fTTl ) are the same and read 



«) = 


(u 2 ) 


{uD = 


(u 2 ) 


[u r Ug) = 


(u 2 ) 


U r U<f,) = 


(u 2 ) 



*2 



2N* 

, n* 2 

2N* 1 
, Cl* 

2N* 2 



■ sin v cos ( 



■ sm( 



(30) 



If the Coriolis number (O is not small (fi* > 1), equations 
( 1291 and (1301 1. though different in details, are of the same or- 
der of magnitude. As the r-approximation, which was used 
to derive these equation, probably has the same accuracy, 
the results are practically the same. We focus in the follow- 
ing discussion on the results of Eq. (130b because the case 
of horizontal original turbulence is probably more adequate 
for stably stratified fluids. 



5 Discussion 

A remarkable feature of Eq. d30b is the finite cross- 
correlation (urU,/,). This means that the turbulence trans- 
ports angular momentum in radius. The cross-correlation is 
positive so that the angular momentum is transported out- 
ward. For the case of differential rotation caused by stellar 
spin-down, this angular momentum flux acts in same direc- 
tion as turbulent viscosity, but it is not viscous by nature. 

What we find here is the known A-effect of differential 
rotation theory (Lebedinskii ll941l Riidiger|T989). Figureg] 
illustrates its origin. The positive correlation (urU^) results 
from the Coriolis force action on the original horizontal 
motion. The fluid particles, which move in the direction 
of global rotation, are deflected outwards in radius by the 
Coriolis force. Fluid particles moving in the retrograde di- 
rection are deflected inwards. The product u r u$ is positive 
in both cases. Note that the influence of the Coriolis force 
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Rotation 




Fig. 4 Illustration of the origin of non-viscous angular 
momentum transport by anisotropic turbulence (see text). 



on the original radial motion produces a negative correla- 
tion (u r u<j,) < (Fig.|4|i. The sense of radial angular mo- 
mentum transport is controlled by the turbulence anisotropy. 
Angular momentum is transported inward if radial mixing 
prevails and outward for predominantly horizontal mixing 
(Rudiger 1989 ). In the present case, the anisotropy of hor- 
izontal type is produced by the stable stratification. It may 
be noticed that we do not find a meridional A-effect. The 
meridional flux of angular momentum is of higher order in 
N*~ 2 compared to the radial flux. This is a consequence of 
small anisotropy in the horizontal plane. Therefore, mixing 
in stellar radiation zones does not produce latitudinal differ- 
ential rotation. 

In a star during its spinning-down, turbulent viscosity 
and A-effect both transport angular momentum outward. If 
the turbulent viscosity in the radial direction is estimated 
as v T k t(u 2 ), the ratio of angular momentum fluxes pro- 
duced by the viscosity and the A-effect can be estimated as 



Qr<f> 



"Iff. I 1 



(31) 



r dn 
q ~ ~0 d7' 

where Eq. ( f30b was used. For the upper radiation zone of the 
Sun, Q 2 /N 2 ~ 1CT 5 (see Fig. 1 in Kitchatinov & Rudiger 
2008). For faster rotating young stars the ratio is larger, but 
still well below unity. Therefore, the angular momentum 
transport by the A-effect is much more efficient compared 
to the effect of the eddy viscosity. 

Radial transport of chemical species is slow because it is 
only due to weak radial mixing and it is further reduced by 
relatively intensive horizontal mixing (Vincent, Michaud & 
Meneguzzi 1996 ; Michaud & Zahn 1998 ). This may explain 
why the characteristic times of rotational coupling between 
the core and the convective envelope ( < 100 Myr; Hart- 
mann & Noyes 1987; Denissenkov et al. 120 101 ) in solar-type 
stars is shorter than the characteristic time (~ 1 Gyr; Sku- 
manich [19721 Melendez et al., |20T0l Baumann et al. |20T0| ) 
of Lithium depletion. 

Papaloizou & Pringle ( 119781 ) found that differential ro- 
tation in radiative cores is unstable to r-modes or vortices 
that are global in latitude but of small scale in radius. For 
the estimates below we assume that the turbulent motions 



are global in the horizontal directions and that radial eddy 
diffusion of chemicals is estimated as 



X T « t{u t ) sa t{u ) 



2\(0) 



2N* 



(32) 



where again Eq. ( 1301 ) was used. Lithium in young solar-type 
stars is depleted about ten times in the first billion years 
of their main-sequence life (Melendez et al. 2010). Lithium 
has to be transported over a short distance (~ O.li?) below 
the convection zone in order to be destroyed (see Fig. 1 of 
Rudiger & Pipin 2001 ). This leads to an estimated fractional 
Lithium abundance (relative to the primordial abundance) at 
the age of 1 Gyr 



exp 



10 9 yr x Xi 
10- 2 R 2 



0.1. 



(33) 



Taking the value of R m 5 x 10 10 cm for the radius of the 
base of (solar) convection zone, we get from Eq. d33l > 



Xt 



2 x 10 3 cm 2 s" 1 . 



(34) 



Equations (l32l and ( l34t lead to an estimate of the character- 
istic eddy turnover time of 

r w 10 7 s (35) 

for a young star rotating ten times faster than the Sun. With 
this value of t, an estimate for the characteristic time of 
core-envelope rotational coupling due to the A-effect can 
be obtained, 



To 



R 2 n 

Qr<f> 



T N 



60 Myr, 



(36) 



where again Eq. (|30T ) for the non-diffusive flux of angular 
momentum = (urU^) was used together with the esti- 
mate N « 10~ 3 s _1 for radiation zones. 

The estimates suggest the following scenario for angular 
momentum and chemical species transport in stellar radia- 
tion zones. Angular momentum loss by a stellar wind in- 
duces differential rotation with an angular velocity decreas- 
ing outwards. The differential rotation is hydrodynamically 
unstable and this instability produces turbulent mixing. The 
stable stratification makes the turbulence highly anisotropic 
with very small radial velocities. The eddy diffusion in ra- 
dius lead to a slow decrease of the Lithium abundance with 
time. Angular momentum is transported not by the eddy vis- 
cosity but by the A-effect due to the correlation of azimuthal 
and radial turbulent motions induced by the Coriolis force. 
The non-viscous outward transport provides relatively fast 
rotational coupling between the radiative core and the con- 
vective envelope. 

Other possibilities for rotational coupling discussed in 
the literature include the angular momentum transport by 
gravity waves (Charbonnel & Talon 2005 ) and by mag- 
netic stress (Charbonneau & MacGregor ll993t Rudiger & 
Kitchatinov 119961 Denissenkov 12010] ). Further studies are 
necessary to decide which of the possibilities are most rele- 
vant. 
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